Synchronization in complex networks with age ordering 
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The propensity for synchronization is studied in a complex network of asymmetrically coupled 
units, where the asymmetry in a given link is determined by the relative age of the involved nodes. 
In growing scale-free networks synchronization is enhanced when couplings from older to younger 
nodes are dominant. We describe the requirements for such an effect in a more general context and 
compare with the situations in non growing random networks with and without a degree ordering. 
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From the brain over the Internet to human society, 
complex networks are the prominent candidates to de- 
scribe sophisticated collaborative dynamics in many ar- 
eas Q. Complex networks are collections of dynamical 
nodes connected by a wiring of edges exhibiting complex 
topological properties. Of particular interest are the so 
called small-world (SW) and scale free (SF) wirings. SW 
are intermediate wirings between regular lattices (RL) 
and random networks (RN) They are characterized 
by a path-length i between any two nodes much shorter 
than a RL, but yet a clustering structure much higher 
than a RN. The basic property of SW is the logarithmic 
scaling of £ with the network size N {£ oc logiV), in con- 
trast to the linear scaling of RL. A specific example of SW 
is the SF configuration, in which the degree k (the num- 
ber of edges in a node) follows a power law distribution 
p{k) . A SF network can be grown by adding suc- 

cessively new nodes to the network and connecting them 
with the already existing ones by a preferential attach- 
ment rule 01 ■ S W and SF properties were found in many 
networks in nature, such as the World Wide Web, the 
internet connection wiring, power grids and transporta- 
tion networks, neural networks, metabolic and protein 
networks. 

Recently, the dynamics of complex networks has been 
extensively investigated with regard to collective (syn- 
chronized) behaviors 0, with special emphasis on the 
interplay between complexity in the overall topologyand 
local dynamical properties of the coupled units US Ql ■ 
Most of the previous works assumed the local units to be 
symmetrically connected with uniform undirected cou- 
pling strengths (unweighted links). This simplification, 
however, does not match satisfactorily the peculiarities 
of many real networks. In ecological systems, for in- 
stance, the non uniform weight in prey-predator inter- 
actions plays a crucial role in determining the food web 
dynamics • Similarly, the interaction between individ- 
uals in social networks @ is never symmetric, rather it 
depends upon several social factors, such as age, social 
class or influence, personal leadership or charisma. 



In this Letter, we analyze networks of asymmetrically 
coupled dynamical units. By explicitly relating the asym- 
metry in the connections to an age order among different 
nodes, we will give evidence that age ordered networks 
provide a better propensity for synchronization (PFS). 
In particular we will show that the three main ingre- 
dient maximizing PFS are i) heterogeneity in the net- 
work topology allowing for the existence of nodes with 
very large degrees (hubs) together with nodes with very 
small degrees (non hubs), ii) asymmetry in the connec- 
tions forcing a preferential coupling direction from hubs 
to non hubs, and Hi) a structure of connected hubs in the 
network. 

We here adopt the idea that the direction of an edge 
can be determined by an age ordering between the con- 
nected nodes. For growing networks (such as SF) the age 
order is naturally related to the appearance order of the 
node during the growing process. We consider a network 
of N linearly coupled identical systems. The equation of 
motion reads 



(1) 



where x = F(x) governs the local dynamics of the 
vector field x^ in each node, H[x] is a linear vecto- 
rial function, and a is the coupling strength. Gij is a 
zero row-sum coupling matrix with off diagonal entries 



Gij — - 
and 0, 
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where A is the adjacency matrix, 

: for i > j {i < j). Ki is the 
set of ki neighbors of the i**^ node, and the parameter 
— 1<6< 1 governs the coupling asymmetry in the net- 
work. Precisely, = yields a symmetric coupling, while 
the limit 6'^— 1(6*— >+l) gives a unidirectional cou- 
pling where the old (young) nodes drive the young (old) 
ones. Asymmetric coupling was recently a,lso established 
for non identical space extended fields 10], where asym- 
metry consisted in forcing preferentially the dynamical 
regime of a field into the other. 

The network PFS can be inspected by linear stability 
of the synchronous state (x^ — Xs,Vi). By diagonaliz- 
ing the variational equation, one obtains N blocks of the 



2 



form (i = JF(xs)(^i + (TAiH[^i], that only differ by the 
eigenvalues of the coupling matrix Ai (here J is the Ja- 
cobian operator). The behavior of the largest Lyapunov 
exponent associated with v — aXi (also called master 
stability function [ij) fully accounts for the linear sta- 
bility of the synchronization manifold. Namely, the syn- 
chronous state (associated to Ai =0), is stable if all the 
remaining blocks, associated to A^ (« > 2), have negative 
Lyapunov exponents. 

For a generic 6, our coupling matrix is asymmetric, and 
therefore its spectrum is contained in the complex plane 
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jX], I = 2, . . . , N). Moreover, since 



all elements of G are real, non real eigenvalues appear 
in pairs of complex conjugates. In the following we will 
order the eigenvalues of G for increasing real parts. Ger- 
schgorin's circle theorem asserts that G"s spectrum 
is fully contained within the union of circles (Ci) having 
as centers the diagonal elements of G {di), and as radii 
the sums of the absolute values of the other elements in 
the corresponding rows ({A/} C U^Gi di,J2j^z I Gij \ )■ 

By construction, the diagonal elements of G are nor- 
malized to 1 in all possible cases. It is crucial to em- 
phasize the physical and mathematical relevance of this 
choice. Physically, this normalization prevents the cou- 
pling term from being arbitrarily large (or arbitrarily 
small) for all possible network topologies and sizes, thus 
making it a meaningful realization of what happens in 
many real world situations (such as neuronal networks) 
where the local influence of the environment on the dy- 
namics does not scale with the number of connections. 
Mathematically, since G is a zero row-sum matrix (and 
furthermore di — '^j^i \ Gij \ because all non zero off di- 
agonal elements are negative), this warrants in all cases 
that G's spectrum is fully contained within the unit circle 
centered at 1 on the real axis (| A/ — 1 |< 1, VZ), giving the 
following inequalities: < A2 < . . . < A^y < 2, and ii) 
I A; |< 1, VL This latter property is essential to provide 
a consistent and unique mathematical framework within 
which one can formally assess the relative merit of one 
topology against another for optimal PFS independent of 
the network size or the local dynamics. 

Let TZ be the bounded region in the complex plane 
where the master stability function provides a negative 
Lyapunov exponent. The stability condition for the syn- 
chronous state is that the set {(7X1,1 = 2, . . . , N} be en- 
tirely contained in TZ for a given a. The best PFS is then 
assured when both the ratio and AI = max/{| X] |} 
are simultaneously made as small as possible. 

With the help of such stipulations, we start with an- 
alyzing the effects of heterogeneity in the node degree 
distribution. This is done by comparing the PFS of a 
class of SF networks with different degree distributions 
with a highly homogeneous RN. The used class of SF 
networks is obtained by a generalization of the preferen- 
tial attachment growing procedure Namely, starting 
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FIG. 1: A5v and (a), 4^(b), and M (c) vs. 6 for SF (m = 5 

and B = Q, solid lines) and RN (dashed lines). All quantities 
and all other network parameters are specified in the text. 
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from TO + 1 all to all connected nodes, at each step a 
new node is added with to links, connecting to old nodes 
with probability pi = ^ '^lk-+B) being the degree of 
the i*'' node, and B a tunable real parameter, represent- 
ing the initial attractiveness of each node). The 7 expo- 
nent of the power law scaling in the degree distribution 
p{k) ~ k~'^iB,ra) -g ^YiQ-a given by 7(B, to) = 3 -h ^ in the 
thermodynamic {N 00) limit 13]. While the average 
degree is by construction (fc) = 2to (thus independent 
of B), the heterogeneity of the degree distribution can 
be strongly modified by B. This induces convergence of 
higher order moments of p{k), in contrast with the case 
5 = that recovers the original preferential attachment 
rule ^J. 

For comparison, a highly homogeneous Erdos-Renyi 
RN is considered, with connection probability P = 
(giving same average degree {k) — 2m), with an 
arbitrary initial age ordering. Fig. ^) shows A^ and 
A2 vs. for SF with m — 5 and B — (solid line) and 
for the chosen RN (dashed line). All calculations refer 
to ensemble averages over 24 different realizations of net- 
works with 500 nodes. For RN, the curve A'^(6') [A^(6')] 
displays a minimum [a maximum] for 9^0, showing 
that asymmetry here deteriorates the network PFS. At 
variance, for SF the difference between A^y and Aj con- 
tinuously shrinks, as 6 decreases. As a consequence. Fig. 
^p) reports the behavior of the eigenratio -r^, making it 
clear that while the best PFS in RN is obtained for the 
symmetric case {9 = 0), SF shows better (worse) PFS 
in the asymmetric case for 9 ^ —1 {9 ^ 1). As for the 
imaginary part of the spectra. Fig. ^) reports M vs. 
9, indicating only very small differences between SF and 
RN in the whole range of the asymmetry parameter, and 
highlighting that the contribution to PFS of the imagi- 
nary part of the spectra does not depend significantly on 
the specific network structure. These findings have been 
consistently observed in subsequent results. 

The second step of our study is the investigation of 
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FIG. 2: 4^ vs. 9 for (a) SF with B = and m = 2 (circles), 

5 (squares), 10 (diamonds); (b) SF with m = 5 and B = 
(circles), 5 (squares), and 10 (diamonds); (c) RN with arbi- 
trary age order (circles), RN with age depending on degree 
(squares), SF with m — 5 and B = without connection 
between nodes 1 and 5 (diamonds), and SF with m — 5 and 
5 = (triangles). 



PFS in SF at different m and B values. In Fig. 13 (a), 
PFS is compared for different m for B = 0. Since m de- 
termines the average degree, as m increases, the average 
connectivity increases and synchronizability is naturally 
enhanced. The most important point is that the mono- 
tonically decreasing behavior of ^ with 9 persists for 
all m, indicating that synchronization is always enhanced 
when 6 becomes smaller. In Fig. 121(b), PFS is compared 
for m = 5 and for various values of B, altering the ex- 
ponent of the degree distribution. For all B values, the 
same enhancement for negative 9 was observed. There- 
fore, one can conclude that in aged growing networks 
PFS depends on the average degree, but the asymmetry 
enhanced synchronization phenomenon does not depend 
on heterogeneity. 

This leads us to discuss the main point of our study, 
concerning the determination of the essential topological 
ingredients enhancing PFS in weighted (aged) networks. 
The first ingredient is that the weighting must induce a 
dominant interaction from hub to non huh nodes. This 
can be easily understood by a simple example: the case 
of a star network consisting of a single large hub (the cen- 
ter of the star) and several non-hub nodes connected to 
the hub. When the dominant coupling direction is from 
the non-hub nodes to the hub node, synchronization is 
impossible because the hub receives a set of independent 
inputs from the different non-hub nodes. In the reverse 
case (when the center drives the periphery of the star) 
synchronization can be easily achieved. The very same 
mechanism occurs in our SF case. Indeed, for positive 
(negative) 9 values, the dominant coupling direction is 
from younger (older) to older (younger) nodes. Now, in 
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FIG. 3: Distribution of G's eigenvalues in the complex plane 
for SF with m = 2, B = and 61 = -0.8 (a), 6* = (b), 
and 9 = 0.8 (c). (d) Master stability function in the complex 
plane for coupled Rossler oscillators. 



SF the minimal degree of a node is by construction m 
and older nodes are more likely to display larger degrees 
than younger ones, so that a negative 9 here induces a 
dominant coupling direction from hubs to non-hub nodes. 

The second ingredient is that the network contains a 
structure of connected hubs influencing the other nodes. 
In our SF case, the normalization in the off diagonal ele- 
ments of G 15] assures that hubs receive an input from a 
connected node scaling with the inverse of their degree, 
and therefore the structure of hubs is connected always 
with the rest of the network in a way that is independent 
on the network size. 

In order to make evident the validity of these claims, 
we have performed a careful analysis on enhanced PFS 
over a series of ad- hoc modified networks. The results are 
summarized in Fig. 121(c). First, we have reordered the 
node age in RN according to each node degree. The re- 
sulting (9) (curve with squares) shows now a minimum 
for an asymmetric configuration (9 « —0.5), in contrast 
to the case with arbitrary aging (curve with circles). This 
confirms the need of a dominant interaction from hubs 
to non hubs for improving PFS, also for highly homoge- 
neous networks. As for the second ingredient, starting 
from a SF with m = 5 and B = (curve with triangles), 
we artificially disconnected the initially existing link be- 
tween the first and fifth network nodes. These are indeed 
the two hubs with highest degree in the SF configuration. 
The result is shown in the curve with diamonds, where 
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FIG. 4: {D) (see text for definition) vs. a for RN (a) and for 
SF (m = 2, _B = 0, b). In both graphs the curves with circles, 
squares, diamonds, up triangles, and left triangles refer to 
0= - 0.8, - 0.4, 0, 0.4, 0.8, respectively. 



one sees that such a small perturbation (the difference 
in the two networks is limited to only a link) is already 
sufficient to substantially weaken the PFS. The situation 
remains however better than the two RN cases, indicat- 
ing that the structure of growing aged network inherently 
enhances synchronization. 

Finally, we prove the validity of our arguments by ap- 
plying them to networks of coupled chaotic Rossler os- 
cillators 1^^. The dynamics is ruled by Eq. ^ with 
x= (a;,y,z),F(x) = [-y-z, x + 0.165y, 0.2-hz(x- 10)], 
and H[x] = x. The master stability function is depicted 
in Fig. 13 (d) in the complex plane ly = a + j(3. The bold 
solid line (denoting a zero Lyapunov exponent) is the 
boundary between stability (7^) and instability regions 
for the synchronization manifold. The whole network 
synchronizes when all the eigenvalues of G (multiplied 
by a) locate inside TZ. Figs. 01 (a,b,c) report the location 
of the G spectrum for SF for m = 2, B = 0, TV = 500 
and with = —0.8, 0, and 0.8, respectively. While for 
9 = Q the spectrum is real, comparison of (a) and (c) 
shows that the spectra at negative 9 values are much less 
dispersed in the complex plane, thus increasing the range 
of a values for which synchronization can be achieved in 
the network. 

The appearance of the synchronous state can be mon- 
itored by looking at the vanishing of the time aver- 
age (over a window T) synchronization error {D) — 

rp(^^_i-^ Sj>i St^^ W'^J ~ ^ilM^'- the present case, we 
adopt as vector norm ||x|| = |a;| -f- \y\ + \z\. Fig. ^re- 
ports {D) vs. 9 for RN with arbitrary age (a) and for 
SF with m — 2, B = (b). The curves with circles, 
squares, diamonds, up triangles, and left triangles refer 
to 6* = - 0.8, - 0.4, 0, 0.4, 0.8, respectively. While 
in the RN case the range for synchronization is substan- 
tially independent on 9 [reflecting the behavior of 



in the dashed line of Fig. ^ (b)] , the case with SF [to 
be compared with the curve with circles in Fig. [21 (a)] 
confirms that synchronization is strongly affected by the 
asymmetry. In particular, negative (positive) values of 
9 have the effects of increasing (decreasing) the range 
of coupling strengths over which synchronization occurs 
with respect to the symmetric case 9 — 0. 

In conclusion we have demonstrated that PFS is en- 
hanced in aged networks of asymmetrically coupled units. 
In growing SF such enhancement is particularly evident 
when the dominant coupling direction is from older to 
younger nodes. Our study allows to individuate the main 
network topological ingredients at the basis of synchro- 
nization enhancement. A key aspect of social organiza- 
tions is the dynamics of information exchange. Our ap- 
proach may provide new insights in the study of collective 
communication or coordination in distributed social net- 
works, as well as useful hints for understanding the for- 
mation of social collective behaviors (leading opinions, 
rumors, political orientations, dominant tastes, habits, 
fashion). Work partly supported by MIUR-FIRB project 
n. RBNE01CW3M-001. 
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